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Abstract 

Quantum tunneling in Reissner-Nordstrom geometry is studied and the tunneling 
rate is determined. A possible scenario for cosmic inflation, followed by reheating 
phases and subsequent radiation-domination expansion, is proposed. 



In 1971, Hawking suggested [T] that there may be a very large number of gravitationally 
collapsed charged objects of very low masses, formed as a result of fluctuations in the 
early Universe. A mass of 10 kg of these objects could be accumulated at the centre 
of a star like the Sun. The masses of these collapsed objects are from 10~ 8 kg and above 
and their charges are up to ± 30 electron units [T|. 

Tracing the evolution of such objects, we propose a mechanism that accounts for the 
cosmic inflation, takes us into a period of reheating phases, and, finally, into the ex- 
pansion of a radiation-dominated Universe. In a nut-shell, the inflation mechanism is 
based on the accumulative effects of Coulomb repulsion at very short range, initially 
completely "cocooned" by Reissner-Nordstrom gravitational effects and subsequently 
unleashed by quantum tunneling. 

Consider the Reissner-Nordstrom geometry [2, 3j in Boyer-Lindquist coordinates 

ds 2 = -4 dt 2 + ^dr 2 + r 2 d6 2 + r 2 S in 2 9 d<\> 2 . (1) 

where: A = r 2 — 2Mr + Q 2 , M is the mass of the centre, and Q — the charge of the 
centre. We will be interested in the case of a naked singularity only, namely: \Q\ > M. 
The radial motion of an ultra-relativistic test particle of mass m and charge q in 
Reissner-Nordstrom geometry can be modeled by an effective one- dimensional motion 
of a particle in non-relativistic mechanics with the following equation of motion [3 El [7] 
(see also [S] for Schwarzschild geometry) : 

r 2 e 2 - 1 

~ + U(r)= € ^, (2) 

where 

C /(r)= 1 (l-4)^-(l-^^ e ) M .-4+ 6 (3) 
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is the effective non-relativistic one-dimensional potential per unit mass, E = (e 2 — l)/2 
is the specific energy of the effective one-dimensional motion, and e = kT/m + 1 is 
the specific energy of the three-dimensional relativistic motion. In equation ([3]), the 
constant a = —Q 2 (l — q 2 /m 2 ) /2 is positive in view of the very high charge-to-mass ratio 
q/m for all charged elementary particles and the parameter b = — M[l — {qQe) / {mM)\ 
depends on the temperature via e. Motion is allowed only when the kinetic energy 
is real. Equation ^ determines the region (r_ ,r + ) within which classical motion is 
impossible. The turning radii are given by El E] : 
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There is no inner turning radius r_ for particles of specific charge q/m such that 
sign(Q)g/m < 1. For particles such that sign(Q)q/m < —1, there is neither inner 
turning radius, nor outer turning radius [21 El E]- Such particles will fly unopposed 
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into the centre. Barrier with two turning radii is present only for particles for which 
sign(Q)g/m > 1. We will consider only such particles. Thus the parameter b will be 
taken as positive. 



U(r) 




Figure 1: The effective potential 

U{r) = | ( 1 — % ~ ~ m M £ ) = ~ ^ + r 

There is no classical analogue of this effect: a charged centre being able to capture 
particles of the same charge within the inner turning radius r_ , despite of the Coulomb 
repulsion. We make the following assumption: the pre-inflationary Universe is an ideal 
quantum gas in thermal equilibrium with constant volume densities of the positive and 
the negative charges. Under minute density fluctuation in the volume density of one 
type of charges at some point, the domain of all other like charges within radius r_ 
are trapped gravitationally into a cluster. We will call this domain a daemon (for dark 
electric matter objects, as introduced by [9], and in line with our previous work [51 16] 
on Reissner-Nordstrom expansion). There will be no charges of this type between r_ 
and r+, while charges of this type on the outside of r + would be strongly repelled. As a 
result, the pre-inflationary Universe nucleates into such domains (daemons). Domains 
of different charge can, obviously, overlap: an oppositely charged particle, approaching 
a daemon, will not experience turning radii, will fly into the daemon and freely interact 
with the particles in it. As our aim is to give a qualitative description, we also assume 
that m is the typical mass of an elementary particle, while q is its typical charge. 
It should be noted that when both turning radii are present, they are always real, 
that is, that they are real for any value of e (or any temperature). The discriminant 
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(4Q 2 /M 2 )(l — q 2 /m 2 )(Q 2 /M 2 — 1) of the quadratic expression in e under the root 
must then be negative and, indeed, it always is — for all charged elementary parti- 
cles, q/m S> 1. Also, an arbitrary accumulation of elementary particles of like charge, 
trapped by the Reissner-Nordstrom field, necessarily leads to \Q\ > M. A daemon is, 
therefore, a naked singularity. 

We now turn to the study of quantum tunneling of trapped particles through the clas- 
sically forbidden region between the two turning radii r±. 

The Schrodinger equation of one-dimensional motion along the r-axis in potential ^ 
is: 

d 2 tp A R 2mE 

^ + (^- 7 )^ = -^W, (5) 

where A = (2m/Ti 2 )a = -mQ 2 (l - q 2 /m 2 )/h 2 = const > and B = (2m/Ti 2 )b = 
-2mM[\ - (qQe)/{mM)]/Ti 2 > 0. 

We are not considering the radial part of a three-dimensional Schrodinger equation as 
we no longer have three-dimensional motion, but an effective one-dimensional problem 
(the difference between the two in our setup is in the parameter B, anyway). 
Essin and Griffiths |10j study very thoroughly the potential — 1/x 2 in quantum mechan- 
ics and its pathologies. They analyze the Shrodinger equation 

Sip a 2mE 

When the positive constant a is smaller than 1/4, there are no bound states {E < 0) 
|10j . (Negative a would turn the potential into a repulsive one.) 

The procedure applied in [TUJ [TT] for < a < 1/4, is not suitable for equation ([5j, due 
to the presence of the 1/r term. Near the origin (r — ► 0), the term proportional to 1/r 
plays little role and asymptotically equation ^ is the same as equation However, 
in view of the parameters involved, q/m S> 1, thus A 1/4, and such situation cannot 
be achieved. 

To address the issue of bound states for equ ation (|5| for the case A > 1/4, we will 
follow the steps of [10]. We introduce k 2 = \J —2mE/Ti. Using the Frobenius method, 
we search for a solution in the form of power series: 

oo 

V(r)=7^5>r\ (7) 

8=0 

where o« = const (with do 7^ 0) and v is also a constant. 
Substituting ([7| into equation (|5]), gives: 

00 

+ v){% + v - X) + A}^- 2 - aiBr 1 - 1 - k 2 ay} = . (8) 

i=0 

The coefficients in the different monomials in r must therefore vanish. 
Setting the coefficient of the 1/r 2 term to zero yields v(y — 1) + A = or 



4 



Here A> 1/4 and thus v is not real. 

Setting the coefficient of the 1/r term to zero gives: 

a\[v{v + 1) + A] - a B = or a\ = — a (10) 

Setting all other coefficients to zero leads to the recursion relation 

a n +2[{v + n + l)(v + n + 2) + A] - a n+1 B - k 2 a n = (11) 
from which a n can be determined in terms of ao for all n > 2 . 

As in the case of [TO] , near the origin (r — > 0), the leading term in the solution is 
a§r v = dQy/r exp[±ig In r] , where g = \J A — 1/4 is real. The solution near the origin 
is real, finite (so that tp — * when r — ► 0) and square-integrable — it is the same as the 
one presented in [10] for equation Q: 

M r) = J 2S ' mH7r9) ^K ig (kr) , (12) 

V 

where Ki g is the modified Bessel function of order ig. 

The allowed energies are not quantized. The problem with this solution is that there is 
no ground state. This means that the particle will cascade down with the release of an 
unlimited amount of energy. 

Regarding scattering states (E > 0), the general solution to the Schrodinger equation 
([5]) near the origin is, again, as the one given in [10] for equation Q: 

Mr) = Vr[FH%\r) + GH$\r)] , (13) 

where k = \?2mE/h , H^' 2 \r) are Hankel functions, and F and G are constants. The 
pathology of this solution is in the fact that the boundary condition at r = imposes 
no constraint on the reflection coefficient and does not determine the amplitude of the 
outgoing wave [TU] . 

Essin and Griffiths [10] propose a renormalization procedure in which the potential 
U(x) = —a/x 2 is replaced by the potential 

TT , . I 00 , for x < a, ,, 

Ujx) = < , 2 r ( 14 

v y [ -a/ or , for x > a . y ' 

The regularized potential has a non-problematic spectrum of discrete bound states and 
a continuum of scattering states [ID]. Upon taking the limit a — > + , all pathologies of 
the — 1/x 2 potential resurface. However, if g also tends to zero, together with a, then 
all excited bound states are "squeezed out" into the continuum of scattering states and 
one single bound ground state with undetermined energy survives [10] . As mentioned 
earlier, in view of the parameters of our daemon model, situation in which g — > is not 
achievable. 

Many other approaches have been considered for resolving the pathologies of the — 1 jx 2 
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potential, in particular, for maintaining a well-defined vacuum. We will only mention 
two — that of Essin and Griffiths [10], in which the non-Hermitian Hamiltonian is 
made self-adjoint by the restriction of its domain, and the one of de Alfaro et al. \FQ , 
in which a different combination of conserved charges was chosen as Hamiltonian. The 
former introduces a free parameter with the dimension of length, thus breaking the 
scale invariance. The latter leads to breakdown of time-translational invariance. While 
a single particle in an —1/x 2 potential does indeed exhibit pathological properties, we 
believe that a proper interpretation can rid us of the pathologies. The situation is very 
similar to the case of a free particle in quantum mechanics. There is no such thing as a 
free particle with definite energy [T3] — the wave function does not represent a physically 
realizable state as it is not normalizable. The physical interpretation of free particles 
comes in terms of wave packets, where interference of waves of different particles leads to 
localization and localization leads, in turn, to normalizability. In the case of the —1/x 2 
potential, the pathologies can be resolved with a likewise many-particle interpretation — 
when packets or, rather, ensembles of particles are trapped by the potential. Consider a 
particle cascading down towards unlimited negative energies. It releases huge amounts 
of energy which, through particle interactions, excite the other particles of the ensemble. 
The excited particles, in turn, stop the cascade of the original particle by exciting it 
with the deposition of positive energy. In other words, the energy exchange between 
different particles does not allow any particle to "shoot down" towards hugely negative 
energies and, in result, the particles "bubble up" at bound states of finite energy. 
In any case, we are interested not in the bound states, but in the continuum of scattering 
states {E > 0). Using the Wentzel-Kramers-Brillouin (WKB) approximation method 
(see [13J, for example), we will determine the transmission coefficient for tunneling 
through the classically forbidden region between the two turning radii r±. The picture 
is very similar to the Gamow theory of alpha-decay (see [13] again) . 
The Schrodinger equation Q5J) can be re-written as: 



where p(r) = ^2m[E — U{r)\ is the classical momentum of a particle with energy E 
moving in potential U{r) (with E > U(r), so that p{r) is real). For tunneling through 
a potential barrier (namely, across the classically forbidden region between the two 
turning radii r±), the WKB-approximated wave function is given by: 



where D = const and \p(r)\ = ^/2m[U (r) — E\. 

The amplitude of the transmitted wave, relative to the amplitude of the incident wave, 
is diminished by the factor e 27 , where 



d 2 i> 



(15) 



dr 2 




(16) 




(17) 
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The tunneling probability P is proportional to the Gamow factor e 27 |13j . 



To solve the integral in (17), we first change the upper limit of integration from r + to 
r + — r, where r is a small positive parameter. Next, we introduce a new variable x, 
using one of Euler's substitutions: 

\/—Er 2 + br — 



(r 



r+ )x . 



(18) 



Therefore, r = (2£r_ + r + x 2 ) / (E + x 2 ) . (The introduction of r does not allow r to reach 



r + where both sides of (18) would vanish and the transformation will not be reversible.) 
Taking the limit r — > + , allows us to express the original integral as in integral over 
the x-axis from to —00. Having in mind that the integrand is an even function of x, 
we get: 



/ 8m E 2 (r + — r_) 5 



r+ 



x 2 dx 



{x 2 + E) 2 {x 2 + Er^/r + ) ' 



The integration gives: 
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Substituting the turning radii Q into this equation gives: 
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m M 
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(19) 



(20) 



(21) 



With the drop of the temperature (that is, when e starts falling from 00 towards 1), 
the inner turning radius r_ tends to a finite value, while the outer turning radius r+ 
tends to infinity. The width of the forbidden classical region, 5 = r + — r_ , also tends 
to infinity in the limit e — > 1: 




Figure 2: Two perspectives of the three-dimensional plot of width of the forbidden classical 
region, 8 — r+ — r_ , as a function of the daemon charge \Q\ and e (both \Q\ and e diminish 
with time). The dashed curve on the left graph indicates a possible history line obtained by 
following a path in the (\Q\, e)-plane. Initially, there is a drop in the width of the barrier 
followed by a rapid increase to infinity. 
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In the very early Universe, at extremely high temperatures (regime e > 1) 
turning radii are approximated by: 

qQ±m\Q\ 



r± 



and 7 is not temperature-sensitive: 



kT 



the two 



(22) 
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7r qQ 
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(23) 



As the emitted particles have charge with the same sign as that of the daemon, the 
absolute value \Q\ of the total charge of the daemon diminishes. The mass M of the 
daemon diminishes as well with each emission, but \Q\/M ~ const > 1 at all times. 
It is, of course, natural to expect that a particle that has just tunneled through the 
potential barrier of one daemon (and emerged on the "outer" side at r+), would tunnel 
through into the "inner" side of a neighbouring daemon with the same charge. This 
happens of course, but not at the rate at which particles tunnel out. In the case of 
standard a-decay, an a-particle is not emitted every time it "knocks" on the "inside 
wall" of the nucleus. The a-particle "rattles" inside the nucleus and the rate of tunneling 
is given by the Gamow factor, multiplied by the factor v/2r, where v is the particle's 
speed and r is the radius of the nucleus [T3] . The "bigger" the nucleus — the smaller 
the rate of tunneling. In our model, three objects with the same sign of their charges 
must be involved in the process of a-like emission and re-capture — two daemons and a 
particle that tunnels out of one and into the other. The two daemons are very strongly 
repelled due to the Coulomb interaction between them. In result, the distance between 
the daemons is much larger than the inner turning radius r_ of each of them. Thus, 
an emitted particle will have to oscillate between two repelling daemons over distance 
much greater than the inner turning radius r_ and the rate of re-capture will be much 
lower than the rate of decay. We will disregard the effect of re-capture of ejectiles on 
the rate of decrease of \Q\ and M. 




Figure 3: Two perspectives of the three-dimensional plot of the tunneling probability P as 
a function of the daemon charge \Q\ and e. Again, both \Q\ and e (or the temperature T) 
diminish with time. The dashed curve on the left graph indicates a possible history line. 
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We next expand the first term on the right-hand-side of equation (23) up to first order 
over the small parameter m/q. This gives the probability for tunneling P as proportional 
to exp[— {ix /K)m^fm \Q\/\q\] in the very early Universe (regime e > 1) and growing 
exponentially with the decrease of \Q\. Over (dimensionless) time dt, the charge of the 
daemon will decrease by the amount d\Q\ proportional to —Pdt and, therefore, in the 
very early Universe, |Q(t)| — ln(C — t), where C = const and t is dimensionless time. 
This gives P(t) ~ 1/(C — i). In alpha-decay, the daughter nucleus recoils after the 
emission. In view of the analogies between alpha-decay and the current case, we make 
the following assumption. A particle of kinetic energy E inside the daemon, tunnels 
through. Tunneling in itself does not change the energy of the particle (otherwise, it 
would "resurface" at point different from the outer turning radius). The recoil energy 
(needed for conservation of momentum) however, does: the particle's kinetic energy after 
the emission will be E, diminished by the recoil kinetic energy E R . The relativistically 
correct relation between the linear momentum p of particle of rest mass m and the 
kinetic energy E of the particle is given by [14J: p 2 = 2mE + 4E 2 /c 2 and the recoil 
kinetic energy E R is [14]: E R = (m/M)E + 2(1 - m 2 /M 2 )E 2 /M where M is the 
mass of the daughter nucleus. Let us first disregard the relativistic (quadratic in E) 
corrections. Then, if Eq is the total kinetic energy of all particles inside the daemon 
before the first emission and if we denote M/m by n (figuratively, we have n "equivalent" 
ingredients inside the daemon), then after the first emission, the ejectile will have energy 
Ei = Eq/u — Ejp, where E R ^ = E\m/[M — m) = E\j{n — 1). The energy of the first 
ejectile is therefore E\ = Eq{ti — l)/n 2 . At the same time, as a result of the loss of 
m/M of the daemon, the total kinetic energy inside the daemon is decreased from Eq to 
Eq — Eq/u = Eo(n—l)/n. The second ejectile will have l/(n— 1) of this energy, or energy 

(2) 

Eq/h prior to leaving the daemon. After tunneling, its energy E 2 will be E /n - E R \ 

where E$ = E 2 m/(M - 2m) = E 2 /{n - 2). Thus E 2 = E (n - 2)/[n{n - 1)]. The 
energy inside the daemon is decreased from Eq(tl — l)/n to Eq{u — 2)/n. The energy of 
the third ejectile prior to leaving the daemon will be l/(n — 2) of the inside energy, or 
Eq/u. Thus the energy carried away by the third ejectile will be E% = Eq{u — 3) /[n(n — 
2)]. The k th projectile will therefore have energy Ej~ = Eq{u — k)/[n{n — k + 1)] = 
Eo(m/M)[l — \j(Mjm — k + 1)]. We now take a continuum limit and re-write this as 
E(t) = Equi/M — Eq/ '[(m/M) {m/M — k(t) + 1)], where k(t) is the number of particles 
emitted after time t. The charge inside a daemon decreses in time from its initial value 
Qo as Q(t) = Q -k{t)q. Thus k{t) = Q /q - (1/q) ln(C- 1) ~ M/m - (l/q) ln(C - t). 
This gives E(t) ~ E {m/M){l-l/[l + (l/q) In(C-t)}}. The temperature drops at least 
as square root of E(t). The outer turning radius r + (which is inversely proportional to 
the temperature) has an accelerated increase with time. In result, the scale factor of 
the Universe, a(t), which is proportional to r + and, therefore, inversely proportional to 
the temperature grows with time. The second derivative of a(t) is positive. Therefore 
we have inflation. 
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Figure 4- Graph of the temperature T as function of time t. The scale factor of the 
Universe, a(t), is proportional to r+ and inveresly proportional to T. The scale factor 
grows with positive second derivative as T decreases. 

If the relativistic corrections p3], mentioned earlier, were included, than the rapid drop 
in T would be even more pronounced. 

Note that in the regime t > 1, the width r+ — r_ = 2m\Q\/(kT) of the classically 
forbidden region initially even decreases with time (as the drop of the temperature T 
is not, initially, as fast as the drop of the charge \Q\ of the daemon — tunneling is 
practically temperature-independent). This is when huge amounts of particles gush 
out of the daemons. The extremely rapid drop in the temperature that follows leads 
to an extremely rapid growth of r+, together with that of a(t). The "graceful exit" 
of the inflation occurs when the width r + — r_ = 2m\Q\/{kT) of the barrier grows 
large enough so that quantum tunneling is switched off. This happens before daemons 
become fully depleted (bound states inside the daemons should also not be forgotten). 
In other words, when the temperature T drops sufficiently and the second term in 



expression \2l\ for 7, namely (vr /Ti)My/m[eqQ / \mM) - l](e 2 - 1)" 1/2 , takes control, 
a break is put on the tunneling (the lower limit of this term is (iT/Ti)qQ/^/rn when 
e > 1). As the probability for tunneling is brought down very rapidly towards and 
particles are no longer ejected by the daemons, the medium outside the daemons is 
no longer cooled by the tunneling process. Without quantum tunneling, the charges 
of the daemons remain practically constant. However, the temperature of the outside 
fraction of the Unniverse continues to drop after the rapid accelerated expansion as a 
different expansion mechanism has naturally taken over. This is the recently proposed 
Reissner-Nordstrom expansion mechanism [5l|6]: with constant charges of daemons, the 
Universe continues to cool: T ~ i -1 / 2 , and expand: a ~ i 1 / 2 . This is the start of the 
radiation dominated epoch. It is also characterized as the beginning of a supercooling 
phase. At the end of the inflation, the daemons are still much hotter than the outside 
fraction of the Universe. A daemon will now cool not through quantum tunneling, but 
through interaction with the particles of oppositely charged daemons, which, in turn 
interact with the particles outside the original daemon. In view of the low densities, this 
does not happen as fast as the Universe expands. Eventually, the temperature of the 
daemons and the temperature of the "free" fraction of the Universe will equalize and, in 
result, the Universe will have reheated, but not enough to reignite the inflation (as the 
daemon temperature now is lower than the one at the end of the inflation and quantum 
tunneling cannot start). During the reheating, the scale factor a(t) of the Universe does 
not decrease as there is no mechanism to draw particles, blown away by the growth of 
the daemons' outer radii, back towards the daemons: the decrease in the outer turning 
radius r + of a daemon simply means that particles of the outer fraction will penetrate 
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deeper and deeper into the repulsive field of the daemons. The Universe then enters 
into another supercooling phase followed by another reheating. This process is repeated 
until daemons cool down to the temperature of the surrounding fraction and cannot 
re-ignite futher reheatings. Then the temperature drop will simply follow T ~ i -1 / 2 
and the expansion will be at the rate of yi. 




t 

Figure 5: Thermal history of the Universe according to the proposed model. The inflation 
(a(i) > 0) is followed by radiation-dominated epoch characterized by Reissner-Nordstrom 
expansion with series of weaker and weaker reheatings after supercooling phases. Overall, 
during radiation domination: T ~ t~ 1/l2 , a ~ t 1 ^ 2 . 

The particles ejected by the daemons, together with all other outside particles, can be 
viewed as a separate fraction of the Universe, additional to that of the daemons. From 
the viewpoint of an incoming non-daemon particle, such that sign(Q)q/m > —1, the 
field of the daemon is characterized by three regions [6] . The first one is the attractive 
region — from infinity down to the radius r c = M(Q 2 /M 2 — 1)[1 — (q/rn)(Q/M)]'~ 1 , 
where attraction and repulsion interchange. There is no gravitationally attractive region 
for incoming particles such that qQ > Mm. The second region is the repulsive region — 
between r c and the outer turning radius r + . Disregarding tunneling back into daemons 
(as described earlier), the third region — between r+ and r = — can be viewed 
as impenetrable. The Universe can therefore be modelled as a van der Waals gas |6 
in view of the deep analogies between the physical picture behind motion in Reissner- 
Nordstrom field and the classical van der Waals molecular model: atoms are surrounded 
by imaginary hard spheres and the molecular interaction is strongly repulsive at close 
proximity, mildly attractive at intermediate range, and negligible at longer distances. 
The laws of ideal gas should then be corrected to account for the increased pressure, 
due to the additional repulsion, and the decreased available volume, due to the presence 
of the daemons. 

We would like to address now the issue of structure formation. This is, essentially, an 
initial-data problem. While the proposed inflation model can indeed accommodate a 
very wide range of initial data, the need of fine tuning with the quantity and composition 
of matter in the Universe and the nature of primeval inhomogeneities is essential for 
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the study of structure formation. On the other hand, it is important to point out that 
the study of structure formation for the proposed model is analogous to the study of 
density perturbations in a van der Waals gas that has undergone accelerated expansion. 
Replacing the perfect fluid equation of state in cosmology with a van der Waals equation 
of state has been considered by many authors. The van der Waals quintessence scenario 
|15[ [T6] achieves exact accelerated expanding solutions. The van der Waals equation of 
state allows only observed fluids to be taken into account; phase transitions to occur in 
the framework of the same evolution; accelerated and decelerated periods to depend on 
the relative values of the parameters of the state equation with respect to the pressure 
and matter energy density which are functions of time [T7] . Additionally, the van der 
Waals equation of state fits the available astrophysical data with the same accuracy 
|18] as the perfect fluid description which works only for very particular conditions and 
scales and which is just a rough approximation of cosmic epochs capable of describing 
stationary situations where phase transitions (which do occur during the evolution of 
the universe) are not considered. The complicated task of studying structure formation 
in a van der Waals universe is addressed by Capozziello et al. in [IT]. In the redshift 
range for the van der Waals quintessence model, where presumably structure formation 
takes place, the baryons energy density dominates over that of the van der Waals fluid 
that, in this period, is very well approximated by a cosmological constant- like term for 
all values of the model parameters [T7]. In the far past, the van der Waals quintessence 
model, is formally equivalent to the ACDM model with the baryons and the Van der 
Waals dark matter playing the roles of CDM and A, respectively [17 . This result 
suggests that structure formation could evolve in a very similar way. 

It is a pleasure to thank V.G. Gueorguiev and R.I. Ivanov for useful discussions. 
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